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ABSTRACT 

We study spectra generated by Comptonization of soft photons by cold electrons 
radially free-falling onto a black hole. We use a Monte Carlo method involving a fully 
relativistic description of Comptonization in the Kerr space-time. In agreement with 
previous studies, we find that Comptonization on the bulk motion of free fall gives 
rise to power- law spectra with the photon index of F > 3. In contrast to some previous 

studies, we find that these power-law spectra extend only to energies <^ TOqC^. We 
indicate several effects resulting in generic cutoffs of such spectra at several tens of keV, 
regardless of any specific values of physical parameters in the model. This inefficiency of 
producing photons with energies > 100 keV rules out bulk motion Comptonization as a 
main radiative process in soft spectral states of black-hole binaries. The normalization 
of the power law (below the cutoff) with respect to the peak of the blackbody emission 
of the surrounding disc is typically very low, except for models with an overlap between 
the disc and the plasma, in which case the spectra are very soft, r>4. 

Key words: accretion, accretion discs - binaries: general - black hole physics - 
radiation mechanisms: non-thermal - X-rays: stars. 



1 INTRODUCTION 

Bulk motion Comptonization (BMC) was first considered 
in a series of papers by Blandford & Payne (1981a, b) and 
Payne & Blandford (1981). They derived the transfer equa- 
tion in the diffusion approximation for photons repeatedly 
upscattered by cold electrons undergoing converging inflow. 
Refining that approach, Psaltis & Lamb (1997) introduced 
relativistic corrections to that transfer equation, solved by 
Psaltis (2001) in a flat spacetime. Solution of the trans- 
fer equation yields a power-law spectrum with a photon 
spectral index of P = 3 for the velocity profile correspond- 
ing to free fall. Titarchuk & Zannias (1998) and TuroUa, 
Zane & Titarchuk (2002) showed that this property remains 
valid when effects of general relativity are taken into ac- 
count. While the above studies demonstrate the ability of 
this model for producing power-law spectra by scattering 
on the bulk motion, they do not determine how far such a 
spectrum extends to high energies. In particular, the Comp- 
ton recoil and other effects resulting in high energy turnover 
of the spectrum were not taken into account in those papers. 

Notwithstanding the uncertainty about that issue, 
Chakrabarti & Titarchuk (1995) and Ebisawa, Titarchuk & 
Chakrabarti (1996) proposed the BMC to be responsible 
for producing high-energy tails in the soft spectral states 
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of black-hole binaries. They argued that both Keplerian 
and sub-Keplerian (close to free fall) components are always 
present in black-hole accretion flows. In the hard state, the 
latter is hot, and thermal Comptonization dominates. In the 
soft state, the sub-Keplerian flow remains cold and Comp- 
tonization on the bulk motion explains the presence of high- 
energy tails, which are observed to have typically F ~ 2.5- 
3.5. Then, application of their BMC model to observational 
data of several sources was presented by, e.g, Shrader & 
Titarchuk (1998, 1999) and Borozdin et al. (1999). 

On the other hand, e.g., Zdziarski (2000), Zdziarski 
et al. (2001), McConneU et al. (2002) and Zdziarski & 
Gierlihski (2004) pointed out that the high-energy tails ob- 
served in the soft states of black hole binaries extend to 
energies well above ~0.5 MeV (see also Grove et al. 1998; 
Tomsick et al. 1999; Ueda et al. 2002), without showing any 
signature of a high-energy cutoff. In particular, the high- 
energy tail of the black-hole binary Cyg X-1 in the soft state 
shows no cutoff up to 10 MeV (McConneU et al. 2002). 

Another disagreements between theoretical predic- 
tions of the BMC model and observational data, noted 
by Zdziarski (2000), concerns Compton reflection (e.g., 
Magdziarz & Zdziarski 1995) which is strong in the soft 
state of black- hole binaries (e.g., Gierlinski et al. 1999), 
while it should be close to null according to the BMC model. 
Furthermore, the presence of the tail in the soft state was 
claimed to be a unique black hole signature as the BMC 
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process requires the presence of a horizon (e.g., Laurent & 
Titarchuk 1999, hereafter LT99). However, by now, observa- 
tions of weakly-magnetized accreting neutron-star binaries 
in the soft state show high energy tails to be common (e.g., 
Di Salvo et al. 2000, 2001, 2002; D'Amico et al. 2001; laria 
et al. 2001; FarineUi et al. 2005). 

Still, the cutoff energy seems to be a key feature allow- 
ing to test the model against data. However, its determi- 
nation has remained relatively uncertain. The studies deal- 
ing with the position of the high energy cutoff include an 
approximate analytic analysis in Titarchuk, Mastichiadis & 
Kylafis (1997) and Monte Carlo simulations by LT99 and 
Laurent & Titarchuk (2001). The Monte Carlo calculations 
of LT99 have been then compared with the CGRO/OSSE 
soft-state spectra of the black-hole binaries GRO J1655-40 
and GRS 1915-)-105 by Zdziarski (2000) and Zdziarski et al. 
(2001), respectively, showing that the BMC predicted spec- 
tra fall below of the observed ones at energies > 100 keV. 

Titarchuk et al. (1997) have presented a simple formula 
for the cutoff energy derived under the assumption that the 
cutoff occurs at the energy for which energy gains are bal- 
anced by energy losses in the plasma rest frame, but ignor- 
ing light bending as well as photon trapping. We find that 
their formula significantly overestimates the cutoff energy 
even if these effects are not important. LT99 use a fully rel- 
ativistic description of photon propagation and scattering in 
the Schwarzschild metric. However, their model appears to 
oversimplify the treatment of photon transfer in an accreting 
plasma, as we discuss in Section |5| 

In this paper, we study formation of BMC spectra us- 
ing a Monte Carlo method described in Niedzwiecki (2005, 
hereafter N05) in a model with radially free-falling electrons 
in a plasma around a black hole. The plasma Comptonizes 
seed photons from an optically-thick disc. Our model in- 
volves a fully general relativistic (GR) description of photon 
transfer and Compton scatterings in a plasma located close 
to a Kerr black hole. In addition, for illustrative purposes, 
we calculate spectra from a model with photon trajectories 
approximated by straight lines. 



2 THE MODELS 

We consider a black hole accreting matter at a mass accre- 
tion rate, M. The black hole is characterized by its mass, 
M, and angular momentum, J. In our model we use three 
reference frames. Photon trajectories are described in the 
Boyer-Lindquist coordinate system, — {t, R, 9, (j)) (which 
generalizes the Schwarzschild coordinate system for J 7^ 0). 
Physical processes (in particular Compton scattering) are 
described in local rest frames of the matter. Finally, refer- 
ence frames of the locally non-rotating observers (observers 
with constant r and 6 but dragged in the azimuthal di- 
rection; in the Schwarzschild metric they are equivalent to 
static observers, i.e., observers with r,9,<l)—const), intro- 
duced by Bardeen, Press & Teukolsky (1972), are auxiliary 
in our model. We use the following dimensionless parame- 
ters; 
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where Me = 4irGMmp/{aTc) is the Eddington accretion 
rate and Rg — GM/c? is the gravitational radius. The for- 
malism used by us does not work for a = (specifically, the 
turning point in ^-motion is found according to the formal- 
ism described in Chandrasekhar 1983, which requires o ^ 0); 
therefore, spectra for a non-rotating black hole are obtained 
assuming a = 10""^. 

For most of spectra presented in this paper, we assume 
that a Keplerian, optically thick, disc is replaced by a free- 
falling plcisma within rtr. We assume = 20, for which a 
relatively high fraction of disc photons is scattered by the 
plasma. The plasma is assumed to be spherical, neglecting 
its likely flattening. Both these assumptions maximize the 
importance of the BMC for spectral formation. In Section 
18.61 we consider also a model where the disc (surrounded by 
the free-falling plasma) extends down to the event horizon. 

Velocity field of the free-fall, with null angular momen- 
tum, is given by 
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yielding the velocity in the locally non-rotating frame. 
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with the corresponding Lorentz factor, 7 = [l — (Z?*^)^] 

We assume that the rest density of electrons, n, is uni- 
form on surfaces of constant r, and determine the density 
radial profile from the continuity equation. Assuming a pure 
H plasma, where the comoving mass density is p = nnip, and 
integrating the mass conservation equation, we obtain, 

" 1 — iJZ' 

R^aT |ff|i/2|u'-|d6l 

where g is the determinant of the metric, \g\ = E^sin^e. 
Below we also use a dimensionless density parameter, 



n = nRsGT. 



(6) 
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For a = 0, equations @ and ^ yield = -(2/r)^/^ and 
h = {2/r)^^^m/2r, respectively, and for a = 0.998, depar- 
tures from these profiles are negligible for r > 3. 

We assume that the electron temperature, Te, is con- 
stant in the inner flow. For most of spectra presented in this 
paper, we assume kTc = 5 keV, at which the BMC domi- 
nates over thermal effects (e.g., Blandford & Payne 1981a; 
Ebisawa et al. 1996) . 

In all simulations, seed photons come from blackbody 
emission of the Keplerian disc. Except for a model with the 
free- falling plasma overlapping with the disc (Section 13.61 . 
we take into account emission from the area of the disc be- 
tween rtr and rout = 100. The latter value is chosen because 
illumination of the inner cloud by emission from r > 100 is 
negligible and thus an increase of rout above 100 does not 
affect resulting BMC spectra. The thermal emission of the 
disc is modelled as in N05; namely, the point of emission of 
each photon is generated according to the radial distribu- 
tion of the disc emissivity (Page & Thorne 1974) and the 
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photon energy (in the disc local rest frame) is generated 
from the blackbody distribution corresponding to the local 
temperature. Note that the radial emissivity (and hence the 
temperature profile) is derived with the boundary condition 
of the vanishing stress at Vms (where rms is the radius of the 
marginally stable orbit), valid for a disc extending down to 
< Tins - In models with the disc truncated at Vti > rms, devia- 
tions from this radial dependence may occur, e.g., a decrease 
of the energy emitted by the disc due to a lower amount of 
mechanical energy transported outward. On the other hand, 
magnetic stresses may be important near the Vms, modifying 
somewhat the above boundary condition and thus the emis- 
sivity profile (e.g., Merloni & Fabian 2003 and references 
therein). However, quantitative description of these effects 
would require hydrodynamical models of the disc and its 
transition into the free-falling plasma, which is beyond the 
scope of this paper. 

Photons illuminating the region of the free-fall are 
traced through their consecutive scatterings until they are 
captured by the black hole, hit the disc, or escape the flow. 
Compton refiection and reprocessing by the disc are not 
taken into account. The spectra presented here correspond 
to the escaping photons as seen by a distant observer and 
are angle-averaged. Simulation of Compton scattering is per- 
formed, in the plasma rest frame, according to the proce- 
dure described in Gorecki & Wilczewski (1984). In modelling 
Comptonization spectra, we consider two cases, one with 
photon trajectories approximated by straight lines, and the 
other with a fuU GR description of photon motion, hereafter 
referred to as the flat and GR models, respectively. 

Our GR model fully follows the procedure described 
in N05, involving solution of equations of photon motion 
in a curved space-time. The optical depth along a photon 
trajectory in a radially infiowing plasma is given by. 



dr — h 



(7) 



This is derived, similarly as in N05, from the elementary ex- 
pression for the probability of scattering in the plasma rest 
frame, dr = nadl, where dl is the length of photon trajec- 
tory measured in the plasma rest frame. In the above equa- 
tion, is an afflne parameter, a is the Klein-Nishina cross- 
section averaged over the Maxwellian distribution of electron 
velocities [see eq. (8) in Gorecki & Wilczewski (1984)] and 
-Erost is photon energy in the plasma rest frame (the energy 
is affected by both the gravitational shift and Doppler shift 
related to the free-fall velocity). Note that increase of the op- 
tical depth in equation Q is written in terms of the increase 
of Boyer-Lindquist coordinates (which describe motion rel- 
ative to a distant observer) along the photon trajectory; the 
derivatives dr/d^ and dt/d<^ are determined by equations of 
motion, see eq. (8) in N05. 

On the other hand, photon trajectories in the fiat model 
are approximated by straight lines and the optical depth has 
the form obtained from equation for r 1, 



dr = na{Ercst)j (1-/3*^ cos (p,-) dl' , 



(8) 



where (f)r is the angle between the outward radial direction 
and the direction of photon motion in the static frame and 
dl' is the distance measured by a static observer (in the 
fiat model, reference frames of static and distant observers 
are equivalent). Note that in our derivation of dr, basing 



on the probability of interaction in the plasma rest frame, 
the kinematic term, 7(1 — (3'^ cos (j)r) (e.g., Blumenthal & 
Gould 1971; Weaver 1976), results from transformation of 
the length of a photon path between the plasma rest frame 
and the frame of a static observer. This term has to be 
taken into account in Monte Carlo simulations of relativistic 
Comptonization (e.g., Pozdnyakov, Sobol' & Sunyaev 1983; 
Gorecki & Wilczewski 1984). Equation Q can be also de- 
rived from the probability of interaction written in terms of 
quantities defined in the frame of a static observer (denoted 
by an apostrophe), dr = a'dl' , where a' = a'{E')n' , taking 
into account the Lorentz invariance of Ea (e.g., Rybicki & 
Lightman 1979, section 4.9). 

All the remaining assumptions in the flat model are the 
same as in the GR model, in particular, capture of photons 
crossing the event horizon (assumed in the flat model at 
rhor = 2 for a = 0). Obviously, our fiat model is not self- 
consistent. We consider it in order to illustrate the crucial 
impact of space-time curvature on emerging spectra. Also, 
a number of previous studies neglected the curvature of the 
space-time and our flat model allows to directly compare 
these previous results with our Monte Carlo simulations. 

We emphasize that an oversimplifled method of inte- 
grating the optical depth along photon trajectories, neglect- 
ing both the gravitational and kinematic effects, have been 
used in some studies of BMC. See, e.g., eq. (2) in LT99 
for the radial Thomson optical depth, which follows directly 
from dr = axndi?, with n being the rest density and ax giv- 
ing the probability of interaction with electron at rest but 
dR being the change of the Schwarzschild radial coordinate. 

In order to illustrate the signiflcance of the relativistic 
terms in equation (|7J , we consider now the Thomson optical 
depth for a radial trajectory in the Schwarzschild metric. 
For a = 0, To = and null angular momentum of a photon, 
equation Q yields 



dr = n7(l ± /3'')(1 - 2/r-) 
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dr. 



(9) 



The last two terms give the proper length in the 
Schwarzschild geometry, gH^dr, where Qrr = (1 — 2/r)~^, 
and 7(1 ± 13"^) gives the Lorentz transformation from the 
plasma rest frame; specifically, 771 gives the density of elec- 
trons in the frame of a static observer and (jt {H l3^) gives 
the probability of interaction with a beam of electrons mov- 
ing with the velocity 13'^ (e.g.. Weaver 1976). For the free-fall 
velocity field, 13^ = and then equation Q yields 

ndr 

The upper and lower sign in equations Q and HIO^ corre- 
sponds to ingoing (dr/d(" < 0) and outgoing photons, re- 
spectively. (Note that the optical depth remains finite for 
ingoing photons at the event horizon.) The non-relativistic 
(for, in particular, /J"" ^ 1) approach, dr = ndr, yields 
roughly two times higher optical depth than our equation 
lp|l for ingoing photons, while for photons escaping from 
the vicinity of the event horizon, it underestimates r by a 
factor of a few. Equation IIUII implies a significant decrease 
of the probability of escape from the innermost region of a 
free-falling fiow with respect to photon transfer in a static 
plasma with the same rest density. 

Monte Carlo simulations of BMC for both fiat and GR 
(with a — 0) models were presented in LT99. In both cases. 
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Figure 1. Emerging spectra for models with a Keplerian disc outside rtr = 20 and a spherical, free-falling, cold (fcTc = 5 keV) inflow 
inside rtr ■ The unscattered part of the blackbody disc emission and the Comptonization spectra are shown by the thin (blue online) and 
heavy (black) solid curves, respectively. All spectra are angle-averaged and correspond to a black hole with M = IOMq at a distance of 
d = 5 kpc. The accretion rate is m = 2, 4, 8 and 12 from left to right. The top panels show spectra for the flat model, and the bottom 
panels are for the GR model with a non-rotating black hole. The dashed (red online), dotted (orange online) and dot-dashed (green 
online) curves show contribution from photons for which the smallest scattering radius is within r = 2—5, 5—8, and 8-20, respectively. 
The (magenta online) points in (e) show the spectrum for rii = 2 and kTc = 5 keV obtained by LT99 in their GR model. 



the assumptions underlying their simulations are similar 
as in this paper (they assume different spatial distribution 
of the source of soft photons, which, however, affects only 
slightly the BMC spectrum, see Section F^.til . However, both 
their flat and GR models appear to neglect the dependence 
of the probability of interaction on the angle between pho- 
ton momentum and the direction of the bulk motion. As 
we discuss above, this may result in strongly overestimated 
emission from the inner region. In Section [3.51 below, we in- 
dicate this effect as a possible explanation for discrepancies 
between results of our model and LT99 at high accretion 
rates. On the other hand, the GR model developed in LT99 

1/2 

takes into account the GR term, g^r , for the proper length 
in calculation of r (L. Titarchuk, personal communication). 

In order to thoroughly test our method and results, we 
have also developed an independent model for Comptoniza- 
tion in the Schwarzschild metric (Niedzwiecki & Sitarek, in 
preparation), which involves a much simpler description of 
GR effects than our Kerr metric model. We have obtained re- 
sults fully consistent with these presented here for a = W~^, 
in particular those where we encountered some discrepancies 
with previous studies. 



3 RESULTS 

Spectra emerging in models with rtr = 20, kTc = 5 keV, 
M = IOMq and various accretion rates are shown in Fig. 
0|a-d) for the flat model and in Fig.0^e-h) for the GR model 
with a non-rotating black hole. We also show contributions 
from photons scattered in various parts of the inner flow. 

3.1 The cutoff energy in the flat model 

In the flat model, the cutoff energy decreases with increas- 
ing m, starting from ~100 keV at low accretion rates to 
a few tens of keV at rh> 10. As illustrated in Fig. ^a~d), 
this decrease of the cutoff energy is primarily related to the 
diminishing contribution of radiation emerging from inner 
parts of the flow. In the flat model, the innermost parts 
generate radiation extending to highest energies. However, 
the optical depth for photons moving outward in the inner- 
most region is typically a few times higher than for photons 
moving inward, see Section |21 As a result, photons at high 
accretion rates are trapped in the inner flow and carried un- 
der the event horizon. Such trapping has been considered. 
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e.g., by Begelman (1979) and Payne & Blandford (1981), 
who defined the trapping radius, rtrap ~ rh, as the distance 
at which the diffusion velocity equals the flow velocity. 

An increase ol m results then in an increase ol the trap- 
ping radius and depletion of increasing number of photons 
from the high-energy part of the spectrum. Moreover, pho- 
tons that diffusively escape from the region inside rtrap lose 
a signiflcant fraction of their energy in scatterings, compare 
the dashed and the heavy solid curves in Fig. 0^d) . In the 
flat model, the position of the trapping radius is crucial for 
the dependence of the cutoff energy on the accretion rate. 

We note that the formula derived in Titarchuk et 
al. (1997) for BMC in flat space-time, Ecut = [4/m -|- 
(4/3)(i''^/c)^]meC^ [note that the form of the second term in 
their eq. (9) is valid only for a plasma with u'' = — (2/r)^'^^c; 
here we use a more general form which can be derived from 
Appendix D of their paper] , predicts the cutoff energies a few 
times higher than these found in our simulations even when 
setting v'' = 0. The cause for that discrepancy is that their 
analysis ignores photon trapping and is based on the as- 
sumption that the cutoff occurs at photon energies for which 
energy gains are balanced by Compton recoil. We find that 
although the average increase of photon energy in a scatter- 
ing indeed increases with decreasing m (as a lower density 
results in a higher average velocity difference between scat- 
terings) , the above formula strongly overestimates the cutoff 
energy. Even at low accretion rates, m < 4, for which trap- 
ping is not important, the formula is not consistent with 
our results, yielding unrealistic values of i?cut > 500 keV, 
regardless of the value of the radial velocity. 

3.2 The cutoff energy including space-time 
curvature 

The spectra in the Schwarzschild GR model are cut off be- 
low 100 keV for all values of the m, see Fig. Q^e-h). Thus, 
bending of photon trajectories strongly diminishes the con- 
tribution from the innermost region with respect to the flat 
model. This is due to several effects. 

As explained below, basic properties of photon motion 
imply that for a = all photons escaping the flow from r < 5 
(i.e., avoiding crossing the horizon) are strongly redshifted as 
seen by a distant observer. Therefore, the emission from the 
r — 2-5 region contributes negligibly to the overall spectrum 
in the GR model with rh = 2, in spite of the trapping surface 
being very close to the event horizon, whereas this region 
gives a major contribution to the 20-200 keV spectrum in 
the corresponding flat model, see Figs. 0^ a), (e). 

The observed photon energy, i?inf, is related to the en- 
ergy in the plasma rest frame, -Brest, by 

£i„f = (1 - 2/r)i''"7(l - /J'-cos^O^rcst, (11) 

where "!>,. is the angle between the inward radial direction 
and the photon direction in the plasma rest frame. Taking 
into account that the Lorentz factor for P'' = 
reduces with the gravitational redshift, (1 — 2/r)^''^ (i.e., 
their product is unity), we obtain, 

£i„f = [l + cos Brest, (12) 

which implies that only photons emitted inward in the 
plasma rest frame (<!>,. < n/2) are observed as blueshifted. 
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Figure 2. (a) The dashed (black online) and solid (blue online) 
curves show the angle-averaged spectra of photons at infinity that 
scatter at 3.5-Rg in the Schwarzschild and extreme Kerr geometry, 
respectively, off electrons with the free-fall velocity {v^ —0.75c 
in both cases). The incident photons, with the power-law spec- 
trum shown by the dotted (red online) line, move in the inward 
radial direction, (b) The angular distribution of photons in the 
local rest frame scattered within the inner 5-Rg for rh = 4. The 
number of photons per unit solid angle is given as a function of 
the angle of the photon direction before scattering with respect 
to the radial direction (cos 3>r = 1 corresponds to the inward di- 
rection). The dashed (blue online) and solid (red online) curves 
show photons with the initial energy corresponding to -Einf < 40 
keV and > 40 keV, respectively, in the Schwarzschild GR model. 
The dotted (black) curve is for photons with E^^i > 40 keV in 
the flat model. The curves are normalized to unity at cos #r = 1. 

On the other hand, the half-angle of the cone of avoid- 
ance, as defined in Chandrasekhar (1983), 

tan* = (r/2 - l)^^\r/3 - l)~'(r/6 + l)"'/^ (13) 

corresponds to $r = tt/S [for which the aberration of light 
gives tan<l? = 1/(7|/3^|)] at r = 5.2. Then, all photons escap- 
ing from r < 5.2 must be emitted backward and thus they 
have -EinfZ-Ercst < 1. The above constraint is valid for a — 0. 
For a = 0.998, such an effective redshift of all escaping pho- 
tons is restricted only to emission from r < 2. Note that 
the solution of the GR radiative transfer equation by Pap- 
athanassiou & Psaltis (2001) similarly indicates that, due to 
the space-time curvature, the innermost region contributes 
negligibly to the BMC spectrum. 

Another reason for the inefficiency of BMC is that the 
free-fall velocity is relatively small in the region where most 
of high-energy emission is generated. In the GR model, the 
major contribution to the high-energy part of the spectrum 
comes from photons scattered within r — 5-8, where velocity 
is at most mildly relativistic, ~ 0.5-0. 7c. 

Finally, we point out the importance of an additional 
effect, occurring for photons undergoing scattering in the 
Klein-Nishina regime, which severely reduces emission of 
hard X-rays. Namely, the highest energies are achieved by 
photons scattered toward the center, which then must be 
scattered off the inward direction to escape the flow. How- 
ever, photons are preferentially scattered forward (and sub- 
sequently captured) at high energies. Then, photons which 
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significantly change direction in a scattering, so that they 
can escape, lose much of their energy to the recoil. Further- 
more, the decline of the Klein-Nishina cross section makes 
scattering of high-energy photons less probable. 

The dashed curve in Fig. |5Ja) shows the spectrum at 
infinity from scattering at r = 3.5 of radially incoming pho- 
tons with a power-law spectrum. At <40 keV, the flux of 
escaping radiation is over an order of magnitude lower than 
the illuminating flux due to capture of the scattered photons. 
At >40 keV, the preference of forward scattering combined 
with the Compton recoil result in a sharp spectral cutoff. 

Fig. Elb) shows the angular distribution of photons un- 
dergoing scattering within inner 5i?g at m = 4. The incident 
photons are strongly concentrated along the inward radial 
direction, and this concentration increases with the increas- 
ing photon energy. This radiation anisotropy results in de- 
pletion of hard X-rays from the radiation leaving the inner- 
most region, crucial for the spectrum. Note the related steep- 
ening above ~40 keV of spectra shown by dashed curves in 
Figs. 0^e-f). On the other hand, the concentration of pho- 
tons in the radial direction is less pronounced in the flat 
model and thus the above effect is less important. 

3.3 The BMC luminosity 

In all the models with rtr ~ 20, the luminosity of the Comp- 
tonization component is <0.04 of the total luminosity. This 
fraction drops to <0.01 when the transition occurs at the 
marginally stable orbit, rtr ~ 6. While the BMC compo- 
nent emerging from a spherically symmetric inner cloud is 
roughly isotropic, the flux from a flat disc observed at an 
inclination, i (with respect to the symmetry axis) decreases 
with increasing i as F{i) oc cosi. As a result, the relative 
flux in the BMC component from systems observed edge- 
on is higher, e.g., ~0.1 for cosi = 0.1, than in the average 
spectra (presented in this paper). 

For rtr = 20, ~0.05 of the photons emitted by the disc 
illuminate the free-fall region. The fraction of photons cross- 
ing the free-fall region that are scattered is 40 and 80 per 
cent for rh = 2 and 12, respectively. Also, both the average 
number of scatterings and the energy transferred to photons 
increase with increasing rh. E.g., Comptonization increases 
the energy of photons by a factor of 1.9 and 4.3 for rh = 2 
and 12, respectively. However, an increase of rh also results 
in an increase of the trapping radius and a higher fraction of 
the photons being captured by the black hole, e.g., 20 and 
40 per cent at m = 2 and 12, respectively. Moreover, trap- 
ping affects primarily photons scattered at small r, which 
achieve the highest energies and carry a significant fraction 
of the total energy transferred to photons. Again for rh = 2 
and 12, 20 and 40 per cent of the photons that get captured 
carry 50 and 85 per cent, respectively, of the total power 
in scattered photons. The resulting conclusion is that the 
luminosity in the BMC component is generally lower than 
the (initial) luminosity in the blackbody photons irradiating 
the plasma. 

3.4 Black hole spin 

Both the radial velocity and density profile changes with 
varying a are negligible for the considered spherical free- 
falling plasma. A significant impact of the space-time metric 



on BMC spectra arises only due to properties of photon tra- 
jectories in the vicinity of the event horizon. The number of 
photons escaping from the innermost region increases with 
the increasing a for other parameters unchanged. Further- 
more, the escaping photons have slightly higher blueshifts 
at higher values of a. For example, in the case of a = 0.998, 
the maximum blueshift, Ei^f/Erast ~ 1.5, is obtained for 
photons escaping from 3 < r < 8, while for a = 0, the max- 
imum blueshift increases from 1 at r = 5.2 to 1.3 at r = 8. 
Such differences in photon motion make the steepening of 
the high-energy part of the spectrum generated by scatter- 
ing at r — 3.5 in the extreme Kerr metric less pronounced 
than in the case of the Schwarzschild metric (Fig. |5^) . 

Fig. 13 compares observed spectra at a = and a — 
0.998 for the same accretion rates. In the case of illumination 
by the external disc. Fig. [^a-cl . Comptonization spectra 
depend very weakly on a. The strongest dependence might 
have been expected for low rh, for which photons scattered 
close to the event horizon are not trapped. For rh = 0.5, Fig. 
13 a), the spectra in models with a — 0.998 and a = are 
dominated above 15 keV by emission from r < 5 and 5 < 
r < 8, respectively. However, the related difference between 
these spectra is rather moderate; most noticeably, the cutoffs 
occur at ~30 and ~50 keV for a = and 0.998, respectively. 
Note also that for such low values of rh, the BMC process 
is very inefficient, yielding spectra with r>4. For rh > 4, 
the spectra are dominated by photons for which the last 
scattering takes place at relatively large values of r, where 
the differences between low and high a are small. This makes 
the black hole spin weakly important for the spectrum. 

Strong gravity effects may affect dynamics of Compton 
scattering in the ergosphere of a rapidly rotating black hole 
(Piran & Shaham 1977). We find, however, that such effects 
are not important for the formation of BMC spectra due to 
strongly relativistic radial velocities. Ergospheric emission of 
a cold free-falling plasma is strongly coUimated inward and 
it does not give any significant contribution to the observed 
spectrum. Thus, we confirm a conclusion of N05 that a quasi- 
Keplerian motion is needed to make radiative processes in 
the ergosphere important for spectral formation. 

3.5 Spectral index 

For all spectra, we determine the photon spectral index, F, 
in the 10-30 keV energy range. Below 10 keV, contribution 
from thermal Comptonization with kT^ = 5 keV is impor- 
tant, and significant departures from a power-law occur in 
some models above 30 keV. 

The values obtained by us for m > 4 in the 
Schwarzschild GR model and for rii > 10 in the flat model 
agree with F = 3 of the solution of the transfer equation in 
the diffusion approximation (e.g., Payne & Blandford 1981). 
At low accretion rates, spectra of the flat model are harder, 
r ~ 2.6 for rh < 4, due to the strong contribution from the 
innermost region. This appears to be due to neglecting terms 
(X iv/cf' in the transfer equation of Payne & Blandford 
(1981), which may significantly underestimate the change of 
energy of a photon scattered in a plasma with a relativistic 
velocity (Psahis & Lamb 1997). 

On the other hand, spectra at low accretion rates in the 
GR model are much softer, r>3.5 for m < 2. A steepen- 
ing with decreasing rh also characterizes spectra obtained by 
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Figure 3. The solid and dashed (blue online) curves show spectra for the GR model for a = and a = 0.998, respectively. The thin 
and heavy curves show the emission from the thermal disc and from the BMC, respectively. The rh = 0.5 (a), 2 (b), 4 (c), and rtr = 20, 
M = 1OM0, kTe = 5 keV. The dotted (red online) curves in (c) are for the GR model with rtr = 6, m = 4 and a = 0. The spectra in 
(d) are for the model with the disc extending down to the event horizon described in Section lij.bl with m = 4 and rfj = 20. 



LT99 for their GR model with fcTe = 5 keV (see table 2 in 
LT99). However, their spectra are harder than ours at any 
value of m. At high accretion rates, the difference is moder- 
ate, e.g., r = 2.8 for m = 7 in LT99 vs. F = 3 in our model, 
and we find that it could be accounted for by including the 
kinematic term in the scattering probability (see Section I^J. 
Indeed, we obtain F = 2.8 with that term neglected in equa- 
tion (|7|. The same conclusion is valid for the flat model. 
Table 1 in LT99 gives T = 2 for their flat model with rh = 7 
and fcTo — 5 keV, while our calculations yield F = 2.7 at 
this rii. However, we obtain similarly hard spectra, F — 2.2 
at rh = 4-8, by neglecting the kinematic term in equation 
JSJ, approximately accounting for the difference. 

At low rh, the discrepancies between our models and 
those of LT99 become more significant, and we have found 
no explanation for them. In particular, we obtain F = 3.6 
for rh = 2 whereas LT99 find F = 2.9. Both of those spectra 
are compared on Fig. 1(e). The (unphysical) modification of 
dr described above does not resolve this discrepancy. 

Furthermore, LT99 conclude that the change of the size 
of the region containing free-falling plasma between rtr ~ 6 
and 40 has a weak effect on the spectra. We find an opposite 
property, with the spectra softening significantly when rtr 
drops below 10. For rtr approaching 6 (at a = 0), the spectra 
soften by at least AF = 1 with respect to those shown in 
Fig. as illustrated by the lower dotted curve in Fig. |3^c). 
The reason for that is obvious; as shown in Fig. a major 
part of the high energy emission is generated in the flow 
beyond 6-Rg. Then the Comptonized component becomes 
much weaker and softer if the flow shrinks to rtr = 6. 

3.6 Illumination patterns 

We flnd that spatial distribution of the sources of seed pho- 
tons on the surface of the outer disc does not change the 



shape of the BMC spectrum. By replacing emission from the 
whole outer disc by emission from its inner edge (at rtr), we 
obtain very similar slopes and cutoffs of the BMC compo- 
nent. This is in agreement with previous studies based on the 
solution of the radiative transfer equation, e.g., Titarchuk et 
al. (1997) and Papathanassiou & Psaltis (2001). 

On the other hand, the relative geometry of the disc and 
the free-fall plasma may be crucial for the normalization of 
the BMC emission with respect to the disc blackbody. The 
relatively low normalization found above (as measured by, 
e.g., the ratio of the EFe of the BMC spectrum at the point 
of intersection with the disc blackbody to the EFe at the 
peak of the disc spectrum) is due to the small solid angle 
subtended by the scattering plasma as seen by the outer 
disc. This effect was noted, e.g., by Borozdin et al. (1999), 
who proposed that an overlap between the free-fall plasma 
with the disc, resulting in a stronger illumination, is needed 
to explain certain observations of soft states of black-hole 
binaries by the BMC model. We note that there is indeed 
some observational evidence for the disc extending close to 
the horizon in some black-hole binaries, e.g., from thermal 
disc emission and relativistic Fe profiles (Zhang, Cui & Chen 
1997; Miller et al. 2002). 

Therefore, we consider a model with a Keplerian, 
blackbody-emitting, disc extending down to the minimum 
stable orbit. Below it, the disc material continues to flow 
to the horizon without dissipation (e.g., Cunningham 1975; 
Muchotrzeb-Czerny 1986). We assume that a spherical 
plasma surrounds the disc at radii r < rg, where a half of 
the accreting matter forms the geometrically thin disc and 
the remaining half accretes through the spherical free fall. 
Then, the density of the plasma is determined by equation 
but with a half of the total rh. The disc emits black- 
body photons down to r^is (= 1.23 and 6 for a = 0.998 
and 0, respectively), but below ra, the temperature profile 
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corresponds to a half of the total rh. The disc absorbs inci- 
dent photons also below rms (as it remains optically thick at 
the accretion rates we consider here; Reynolds & Begelman 
1997). 

Fig. Eld) shows the resulting spectra for a model with 
Tff = 20 and the (total) m = 4. The BMC component con- 
tains now 6 and 18 per cent of the total luminosity for a = 
and 0.998, respectively. On the other hand, the spectra of 
this component are significantly softer (by AT ~ 0.5 for 
a — 0) than in the corresponding (with m — 2) model with 
truncated disc. 

For a = 0, >70 per cent of the energy dissipated in the 
disc is released at r > 20. Therefore, a disc truncated at 
r = 20 has a similar luminosity as a disc extending down to 
Tms- However, a significant fraction of soft photons emitted 
within the plasma is scattered in the latter case, which re- 
sults in the fraction of disc photons that are Comptonized 
increasing to 10 per cent (from 3 per cent for a truncated 
disc with m — 4). A related increase of the peak of the BMC 
component at ~1.5 keV (formed mostly by the first scatter- 
ing order) by a factor of ~3 is clearly seen by comparing the 
heavy solid curves in Figs. 3(c) and (d). 

On the other hand, over 30 per cent of the scattered 
photons are absorbed by the part of the disc within the free- 
fall region. Moreover, this absorption affects mostly photons 
scattered in the inner region. In the corresponding model 
with a truncated disc, a significant fraction of photons form- 
ing the high-energy part of the spectrum gets scattered at 
r = 5-8 inward but then escape (after passing the turn- 
ing point in the r-motion), crossing the equatorial plane. In 
contrast, those photons are absorbed in the model with the 
overlapping disc. On the other hand, this effect is much less 
significant for photons at lower energies, typically formed at 
larger distances. Thus, the decreased contribution of high 
energy photons results in a softer spectrum, with F — 4.1 
vs. r — 3.6 in the truncated disc model with the same elec- 
tron density. 

In the case of a = 0.998, the seed photons are very cen- 
trally concentrated, with over a half of the disc luminosity 
produced at r < 5. Thus, they have relatively high energies, 
of several keV. As a result, a power-law shape of the spec- 
trum is achieved only in a narrow energy range of 20-70 
keV. Absorption of photons with highest energies by the disc 
is significant in this photon trajectories are strongly 

bent toward the equatorial plane at small r in the extreme 
Kerr metric. A resulting BMC spectrum is shown by the 
heavy dashed curve in Fig.j^Jd). 

In the overlap model, an increase of m yields a decrease 
of the relative normalization of the BMC component for a = 
0.998, e.g., its contribution to the total luminosity decreases 
from 18 to 10 per cent for m changing from 4 to 24. For 
a = 0, this relative fraction remains approximately constant, 
at ~6 per cent. The fact that the efficiency of the process 
does not increase with increasing m is due to the increasing 
number of photons trapped at high rh fSection l3.3ll . 

The value of is crucial for the normalization of these 
models. A high normalization of the BMC component with 
respect to the disc one may be achieved if a very dense, free- 
falling, plasma extends out to further distances than consid- 
ered above. If an optically-thick free-fall plasma covers the 
disc surface up to > 100, most of disc emission is Comp- 
tonized. We note that the BMC model makes some specific 
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Figure 4. The dotted (green online) and dashed (red online) 
curves show Comptonization spectra from static and free-falling, 
respectively, plasmas at kTc = 50 keV. The heavy solid curve 
shows the spectrum from a free-falling plasma with fcTe = 5 keV. 
The free-fall spectra are for rh = 4 and the static plasma spectrum 
is for the same density distribution. All models have a = 0, rtr = 
20 and M = IOMq. The thin (blue online) solid curve shows the 
blackbody disc emission (the same for all three models). 

predictions for such scenario. Namely, the optical thickness 
of the free-fall plasma in the region within r ~ 20-100 is a 
few times lower than the optical thickness of the inner re- 
gion, r < 20. Then, if the free-falling plasma is sufficiently 
dense to Comptonize most of disc emission from r > 20, al- 
most all of the total emission from the inner region will 
be trapped. Then, any relativistic signals from the inner- 
most region, including relativistic fluorescent lines or high- 
frequency QPOs (assuming their relation to Keplerian fre- 
quencies), would not be observable. Furthermore, the BMC 
photons observed by a distant observer would come from 
scatterings at large r in a plasma with subrelativistic veloc- 
ities. Then, spectrum would be precisely approximated by 
the flat model, i.e., would have F = 3. 

LT99 argue that illumination of the central region can 
significantly increase if the surrounding disc is geometrically 
thick. They approximate this effect by emission from the 
inner edge of the disc at 6i?g. However, emission from the 
surface of the disc is not shown in that paper, and thus the 
actual normalization is uncertain. 

Summarizing this Section, we confirm that the relative 
normalization increase significantly when the source of seed 
photons is within the free-fall region. On the other hand, 
we find that then the BMC spectrum component becomes 
significantly softer than in corresponding models without 
overlap. 

3.7 BMC vs. thermal Comptonization 

Fig. 131 compares a BMC spectrum with spectra for semi- 
relativistic thermal Comptonization. For kTc ^ 10 keV, ther- 
mal Comptonization completely dominates formation of 
spectra and the role of bulk motion is negligible (see also 
Jaroszynski 2001). 

Fig.Elalso compares spectra from thermal Comptoniza- 
tion between static and free-falling plasmas with the same 
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rest density distributions. The major impact of the free fall 
on formation of thermal Comptonization spectra is by trap- 
ping photons in optically thick regions of inflowing plasma, 
making a Wien bump much less pronounced (which effect 
has been studied by Colpi 1988). 



4 DISCUSSION AND CONCLUSIONS 

We have examined formation of spectra emerging due to 
multiple scattering of photons in black-hole accretion flows 
with relativistic radial velocities. In agreement with previ- 
ous studies, we find that this model gives rise to power-law 
spectra with the photon spectral index of F « 3 at high m, 
which soften to F > 3 with decreasing m. The BMC process 
cannot explain spectra with F < 3. At semi-relativistic tem- 
peratures, needed to obtain such hard spectra by thermal 
Comptonization, the BMC process contributes negligibly to 
the spectral formation. 

We point out effects preventing this power-law from ex- 
tending beyond ~100 keV. The recoil effect, dominant at 
lower accretion rates in flat space-time, results in a cut- 
off around 100 keV. At higher accretion rates, trapping of 
photons reduces contribution of photons from innermost re- 
gions. In a curved space-time, bending of photon trajectories 
combined with an increased probability of forward scatter- 
ing at high energies result in cutoffs at a few tens of keV 
even at low accretion rates. 

We do not confirm claims of, e.g., Titarchuk et al. (1997) 
and LT99, that the power- law spectrum extends to ~ rricC^. 
We find that such high energies are achieved only by photons 
scattered toward the center at very small radii, which then 
do not contribute to observed spectra. 

In models with illumination of the free-falling plasma by 
the surrounding disc, the normalization of the BMC spectral 
component with respect to the peak of the disc blackbody 
spectrum is very low, typically < 10~^. If a likely flattening 
of the plasma were taken into account, this relative nor- 
malization would become even lower. Higher normalizations 
may be obtained in models involving a disc extending down 
within the free-falling plasma. In this case, the high nor- 
malization occurs together with soft BMC spectra, F > 4, 
due to absorption of Comptonized photons with highest en- 
ergies by the inner disc. The luminosity of the (very soft) 
BMC component comparable to the disc one can only be 
achieved in models with complete obscuration of the central 
region, implying then no observable relativistic signatures. 

We find that rotation of the black hole has a minor effect 
on the shape BMC spectra. The rotation is only important 
for models with the disc extending close to the event horizon, 
as in this case the fraction of soft seed photons emitted in 
the innermost region increases with an increase of a. 

The amplitudes and spectral indices of high energy tails 
in some observations of black-hole binaries are consistent 
with the BMC model for moderate accretion rates, ?ti<2. 
E.g., according to table 1 in Borozdin et al. (1999), the 1997 
March 24 observation of GRO J 1655-40 is consistent with 
parameters corresponding to the external illumination, and 
observations of 1997 April 24 and 30 are consistent with the 
untruncated disc. However, as discussed in Section Q deter- 
mination of the position of the high-energy cutoff is crucial 
for the confirmation that those spectra may be explained 



by the BMC model. As shown by Tomsick et al. (1999), the 
high-energy tail of this object extends to ~2 MeV, i.e., much 
higher than that possible to obtain in this model. In general, 
the BMC model is unable to explain the high energy tails of 
black-hole binaries in the soft state, which have often been 
measured to extend to energies > 500 keV without any cut- 
off, and often have too high amplitudes and/or F < 3 (see, 
e.g., Zdziarski & Gierliriski 2004). 

A very low normalization of the high-energy tail (with 
a poorly constrained high-energy cutoff) is often seen in 
the ultrasoft (i.e., strongly dominated by a disc blackbody, 
e.g., Gierlinski & Done 2004) states of black-hole binaries, 
see fig. 8 in Zdziarski & Gierliriski (2004). However, those 
tails appear rather hard, e.g., with F 2 in GRS 1915-1-105 
(Zdziarski et al. 2001). 

Finally, we caution that a commonly-used implementa- 
tion of the BMC model, BMC in the X-ray data fitting code 
XSPEC, allows for an arbitrary normalization of the Comp- 
tonization tail, arbitrary value of the spectral index as well 
as it does not include any high-energy cutoff, in contrast to 
the theoretical results. Thus, satisfactory results of fitting 
X-ray data with this model cannot be automatically taken 
as showing a consistency between that model and the data. 
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